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Abstract — Recently, it has been demonstrated experimentally 
that adaptive estimation of a continuously varying optical phase 
provides superior accuracy in the phase estimate compared to 
static estimation. Here, we show that the mean-square error 
in the adaptive phase estimate may be further reduced for the 
stochastic noise process considered by using an optimal Kalman 
filter in the feedback loop. Further, the estimation process can 
be made robust to fluctuations in the underlying parameters of 
the noise process modulating the system phase to be estimated. 
This has been done using a guaranteed cost robust filter. 

I. INTRODUCTION 

Quantum parameter estimation is the problem of esti- 
mating an unknown classical parameter, often an optical 
phase shift, of a quantum system [1]. It is at the heart of 
many fields such as gravitational wave interferometry [2], 
quantum computing [3] and quantum key distribution [4]. 
The fundamental limit to the accuracy of the optical phase 
estimate is set by the Heisenberg's uncertainty principle [5], 
By contrast, the standard quantum limit (SQL) refers to 
the minimum level of quantum noise that can be obtained 
using standard approaches to phase estimation which do 
not involve real-time feedback. The SQL sets an important 
benchmark for the quality of a measurement and provides an 
interesting challenge to devise quantum strategies that can 
beat it. 

Up until very recently, most of the work in quantum phase 
estimation has been on estimating a fixed unknown phase 
shift. It was shown theoretically that adaptive homodyne 
single-shot measurements can yield an estimate with mean- 
square error less than the SQL [6], [7], [8]. This was 
subsequently demonstrated experimentally using very weak 
coherent states [9]. However, a more experimentally relevant 
problem is when the phase varies continuously under the 
influence of an unmeasured classical stochastic noise process 
[10], [11], [12]. 

The first experimental demonstration of adaptive quantum 
phase estimation of a continuously varying phase was pre- 
sented in Ref. [13], where an estimate could be obtained with 
a mean-square error of up to 2.24 ± 0.14 times smaller than 
the SQL. In the adaptive experiment, the system phase to be 
estimated was modulated by a classical stochastic Ornstein- 
Uhlenbeck (OU) noise process. We show here that the mean- 
square error in the estimate can be further reduced for the 
case of an OU noise process by using an optimal Kalman 
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filter in the feedback loop and that the filter used in Ref. 
[13] is only optimal for the case where the noise process 
is a Wiener process and the measurement is assumed to be 
linear. See also Ref. [14]. 

It is physically unreasonable to specify precisely the de- 
sired values of the underlying parameters of the noise process 
modulating the system phase to be estimated. Hence it is 
desired to make the estimation process robust to uncertainty 
in these parameters. A robust quantum parameter estimation 
technique, as applied to atomic magnetometry, was demon- 
strated theoretically in Ref. [15]. Also, a simple approach 
to robust adaptive phase estimation was discussed briefly in 
Ref. [14]. Robust adaptive estimation of the continuously 
varying optical system phase can be made by making the 
Kalman filter in the feedback loop robust to uncertainties 
introduced in the underlying parameters of the noise process. 
In this paper, this is achieved by using the guaranteed cost 
robust filtering approach described in Ref. [16]. 

II. MODEL OF ADAPTIVE PHASE ESTIMATION 
TECHNIQUE FROM REF. [13] 

Fig. Q] shows the model block diagram of the adaptive 
system with filter used in Ref. [13] in the feedback loop. 
The estimator calculates offline the estimate for the system 
phase. 

A. Process 

The governing equation for the adaptive estimation sys- 
tem being considered is the dynamically varying stochastic 
classical phase [13]: 



d(j){t) = -X(f>(t)dt + s/~KdV{t), 



(1) 




Fig. 1. Block Diagram of the adaptive system considered in Ref. [13]. 



where <fi(t) is the system phase to be estimated, dV is a 
Wiener increment, A > is the mean reversion rate and 
k > is the inverse coherence time. 

The process model for the above system is, therefore, 
given by the equation: 



4>(t) = -\<f>(t) + y/Kv{t), 



(2) 



where v := ^Jf is a white noise process with autocorrelation 
function R v (t) = S(r). 

B. Measurement 

Using a linearization approximation, the homodyne pho- 
tocurrent from the adaptive phase estimation system is given 
by [13]: 



I(t)dt = 2\a\[(f>(t) - $(t)]dt + dW(t), 



(3) 



where \a\ is the amplitude of the coherent state with photon 
flux given by J\f := \a\ 2 , cf> is the intermediate phase estimate 
which is also the optical local oscillator phase because of the 
feedback, and dW is Wiener noise arising from the quantum 
vacuum fluctuations. 

The instantaneous estimate 6(t) is given by: 

m 

2\a\- 



6(t) := 0(t) 



(4) 



The measurement model can, therefore, be described by 
the equation: 



o(t) = 4>(t) + J^Mt), 



(5) 



where w := is also a white noise process with R w (t) — 
S(t). 

C. System Model 

Rewriting the equations for the system under considera- 
tion, we get 



Process model: 4> 
Measurement model: 9 



= —\<t> + \/~KV, 



(6) 



where 



E[v(t)v(r)}=QS(t-r), 
E[w(t)w(r)} = R5(t-r), 
E[v(t)w(r)] = 0. 

Since V and W are Wiener processes, both Q and R are 
unity. 

D. Feedback 

1 ) Transfer Function: The feedback filter used for adap- 
tive phase estimation in Ref. [13] has the following form: 

ft 

(7) 



e-(t) = x/ e{T)e x < T - t >dT = x (e(t) 



,-x-f) 



where the value of \ is \opt — 2\a\^/~K, which is optimal in 
the limit A -t 0. 



The transfer function of this filter can, thus, be found to 



be: 



G P ( S ) 



e-QO 

9(8) 



(8) 



2) Error Covariance: We augment the system given by 
(O with the feedback filter $Q as shown in Fig. [T] and 
represent the augmented system by the state-space model: 



= A: 



where 



6- 



and 



(9) 
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From (O and ([8j, we obtain: 

4> = — + ^/kv, 
©- = ~Xopt<d- + Xo P t4> - 



Xopt r 

2M'' 



(10) 

(11) 



Thus, we have: 

-A 

Xopt - 





Xopt 



and B 









Xopt 

2\a\ 



For the continuous-time state-space model (O, the steady- 
state state covariance matrix P5 is obtained by solving the 
Lyapunov equation: 



AP 6 + P S A 
where P5 is the symmetric matrix 



BB T = 0, 



(12) 



P s = £(xx' ) = 
Upon solving (TT2l . we get 



Pi P2 
P2 P3 



P - K 
P2 - 



XoptK 



P:i 



2A(A + xopt) ' 
4-XoptM 2 K + Xopt A 2 + x 2 p t A 



8|a| 2 A(A + x opt ) 
The estimation error can be written as: 

e = (j) - 9_ = [1 - l]x, 

which is mean zero since all of the quantities determining e 
are mean zero. 

The error covariance is then given as: 



a 2 P = E(ee T ) = [1 - l]fi(xx r ) 

1 



= [1 - 1] 
Thus, we obtain 



Pi P2 
P2 P3 



-1 



1 



Pi - 2P 2 



dp 



Xopt (A + 2x pt) _ Vk(A + 4|a|y^) 
8|a| 2 (A + xopt) ~ 4|a|(A + 2|a|V^)' 



(13) 



One can verify that this analytical expression for the error 
covariance agrees with equation (10) of Ref. [13] for the 
optimal case of x- 



III. MODEL OF ADAPTIVE PHASE ESTIMATION 
USING A KALMAN FILTER 

Fig. |2] shows the block diagram of the adaptive system 
with a Kalman filter in the feedback loop. As compared to 
Fig. [U there is a subtle difference in the way the input 9 to 
the feedback filter is generated in this case. 

A. Algebraic Riccati Equation 

Given the system as described by ©, the steady-state 
continuous Kalman filter Riccati equation [17] can be shown 
to be: 

(14) 



-2AP 



4|a| 2 P 2 



K = 0. 



B. Filter Equation 

The Kalman filter equation for the (continuous) system 
under consideration is then given by: 



= -(\ + K)4> + K6, 



where K is the Kalman gain. 

The stabilising solution of ( TT4T > can be found to be: 



P 



4|a 



(15) 



(16) 



The Kalman gain for the system under consideration is 
then given by: 



K 



-A + y 4k|q;| 2 + A 2 . 



(17) 



C. Feedback 



1) Transfer Function: The transfer function of the 
Kalman filter may be obtained from $15[ to be: 



G K {s) 



(s) 



K 



s + X + K 



(18) 



2) Error Covariance: The error covariance for the 
Kalman filter is given by i 161 . rewritten as below: 



' K 



K 



(19) 



One can verify that the error covariance for the Kalman 
filter, obtained using the Lyapunov method used in section 
III-D. 21 earlier, would be the same as in the above equation. 




I = I + 2\a\(f> 



Kalman 




V' 


Filter 






G K (s) 


e 



Fig. 2. Block Diagram of the adaptive system with a Kalman Filter in the 
feedback loop. 



IV. COMPARISON OF KALMAN FILTER WITH 
FILTER USED IN REF. [13] 

The filter used in the feedback loop for the adaptive system 
considered in Ref. [13] was designed to behave optimally 
under the condition that the value of A is zero, so as to 
approximate the detailed analysis in Ref. [10]. Here, we will 
see that a Kalman filter would be optimal in all cases, i.e. for 
all values of A and that the filter used in Ref. [13] coincides 
with the Kalman filter in the limit A — >• 0. 

A. Transfer Function 

Equations (O and ( TT~8b given earlier are the transfer 
functions of the filter used in Ref. [13] and the Kalman filter, 
respectively. 

In the limit A — > 0, we have 



G K {s) 



K 



where 

K - 



Urn (-A + V4«H 2 + A 2 ) = 2\a\V^ = Xo P t- 



Let us consider the parameters for the adaptive phase 
estimation problem considered in Ref. [13]: Map = |c*| 2 = 
1.3499 x 10 6 s" 1 , k ap = 1.5868 x 10 4 rad/s, \ AP = 
6.1451 x 10 4 rad/s. 

We can calculate the following, using the units as indicated 
above: 2\a\ = 2324, = 126, Xo P t = 292824, and K = 
237643. 

Thus, we have: 

292824 

Gp{s) = 



whereas 



G K {s) = 



+ 292824' 
237643 



299094 

The transfer function of the filter used in Ref. [13] is 
similar to, but not the same as, that of the optimal Kalman 
filter for the given set of parameters. 

B. Error Covariance 

Equations (fT3l and ( fl~9T > given earlier are the error covari- 
ances in the estimates from the filter used in Ref. [13] and 
the Kalman filter, respectively. 

Note that in the limit A — > 0, we have 



Co = ov = 



2\a\ 



(20) 



Considering again the parameters for the adaptive phase 
estimation problem considered in Ref. [13], we can evaluate: 



<t 2 p = 0.0495, 



whereas a\ = 0.044. 



Thus, it is clear that the filter used in Ref. [13] is sub- 
optimal as compared to the use of a Kalman filter for the 
given value of A. 

However, in the limit A — > 0, we would have: 

a], = o% = 0.0542. 



Fig. [3] shows the plot of the error covariance against the 
parameter A for the two cases, viz. Kalman filter and the 
filter used in Ref. [13]. Here, we have used the nominal 
experimental values for the other two parameters, viz. \a\ 
and k, in the expressions dT3b and ( fT9b . As can be seen, 
the filter used in Ref. [13] behaves exactly like the optimal 
Kalman filter for lower values of A. However, as the value 
of A rises, the Kalman filter's error covariance improves 
significantly as compared to that of the filter used in Ref. 
[13]. The red vertical line denotes the value of A for the 
adaptive experiment considered in Ref. [13]. 

V. ROBUST FILTER 

In this section, we make our filter robust to uncertainty in 
one of the underlying parameters using the guaranteed cost 
estimation robust filtering approach given in Ref. [16]. 

A. Process and Measurement 

We introduce uncertainty in the parameter A as follows: 

A -)• A - [i\A, 

where A is an uncertain parameter which satisfies |A| < 1 
and < /i < 1 is a parameter which determines the level of 
uncertainty in the model. 

The process and measurement models of (0 take the form: 



Process model: 4> = (—\ + fi\A)(j) + t/kv, 

1 



Measurement model: 9 



2\a\ 



(21) 



B. Riccati Equation and Optimal Error Bound 

As in Ref. [16], the Riccati equation for the guaranteed 
cost filter for the system is: 

/j 2 A 2 



4|a| 2 Q 2 



2AQ 



0. 



(22) 



The stabilising solution of this equation yields an upper 
bound for the robust filter error covariance: 

+ _ -Ae + sj\ 2 e 2 - e 3 K - e 2 fi 2 \ 2 + 4\a\ 2 e 2 K + 4\a\ 2 efi 2 X 2 



e(-e + 4|a| 2 ) 



(23) 
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Fig. 3. Comparison of the error covariance between the Kalman filter and 
the filter used in Ref. [13]. 



The optimum value of e at which the bound Q + is the 
minimum can be found to be: 



+ A + sjti 2 \ 2 - 2\ 2 fi + A 2 + 4|q| 2 k) fiX 



(24) 



C. Filter Equation 

We calculate the robust filter equation for our system to 
be [16]: 



■ = -\$ + (e - 4|a| 2 )Q+0 + 4\a\ 2 Q + 4> + 2\a\Q + w. (25) 



Note that when /i = 0, then e opt = and ( 1221 reduces to 
([T4l with P replaced by Q. That is, the robust filter reduces 
to the standard Kalman filter. 

D. Transfer Function 

Equation d25l > yields the below transfer function for the 
robust filter: 



G R (s) 



(s) 



4|a| 2 Q+ 



6(s) s + A + (4|a| 2 - e)<^ 



(26) 



This transfer function is a first-order low-pass filter with 
gain and corner frequency slightly different from those of 
the filter used in Ref. [13] and the Kalman filter. 

VI. COMPARISON OF ROBUST FILTER WITH 
KALMAN FILTER 

We can compute using the Lyapunov method employed in 
section III-D.2I and for the nominal experimental values of all 
the parameters and /1 = 0.05, the error covariance for our 
robust filter as a function of A. We can similarly compute 
the error covariance as a function of A for the Kalman filter 
(fT~5b where the process and measurement models are for the 
uncertain system given by (1211 . 

We can then obtain a plot of the error covariance as a 
function of A for the robust filter and the Kalman filter on 
the same graph, as shown in Fig. [4] As we can see from this 
figure, there is not a significant performance improvement 
with the robust filter as compared to the Kalman filter for the 
case of 5% uncertainty. The horizontal dashed line indicates 
the optimal error bound for the robust filter. 

We can, similarly, plot such comparison graphs for the 
case of more uncertainty in A. Figs. [5] and [6] show plots for 
H = 0.5 and /.i = 0.8, respectively. As we can see, the robust 
filter performs better than the Kalman filter as A approaches 
1 for all levels of uncertainty in A. As the uncertainty in A 
increases, so does the improvement in performance of the 
robust filter relative to the Kalman filter. 

VII. CONCLUSION 

This paper applies the theory of robust filtering to the 
problem of adaptive homodyne estimation of a continuously 
evolving optical phase shift of a quantum system. The 
immediate further work to this would be to extend the 
theory to include Kalman and robust smoothing rather than 
filtering alone. It remains to illustrate the results in this article 
experimentally. The results herein may as well be extended 
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Fig. 4. Comparison of error covariance as a function of A for /i = 0.05. 




Fig. 5. Comparison of error covariance as a function of A for = 0.5. 
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Fig. 6. Comparison of error covariance as a function of A for fj, = 0.8. 



for the case of squeezed states of light. Also, it would be 
interesting to explore robustness as applied to other types of 
complex noise processes or uncertainties in other parameters 
such as the noise power or photon flux. 
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